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flow  computer  codes  for  rotating  blades  to  form  an  efficient  turbulent 
boundary-layer  computer  code.  For  a  given  potential-flow  solution,  a 
typical  boundary-layer  computation  requires  less  than  10  seconds  computer 
time  on  the  Burroughs  7700  high-speed  computer.  Boundary-layer  predictions 
are  presented  for  several  rotating  blades.  Computed  results  are  shown  to  be 
in  agreement  with  experimental  data  for  a  simple  rotating  body.  For  the 
examples  considered,  displacement  of  the  mid-chord  point  of  a  blade  from  a 
straight  radial  line  is  predicted  to  reduce  the  computed  values  of  local 
skin  friction  coefficient,  with  an  estimated  increase  in  overall  eftfciency 
of  about  one  percentage  point. 
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ABSTRACT 


A  general  formulation  of  the  boundary-layer  computation  scheme  on  the  surface  of  a 
rotating  blade  is  presented.  Momentum-integral  methods,  together  with  the  three-dimensional 
entrainment  equation  for  a  rotating  disk,  are  used  to  calculate  the  three-dimensional  turbulent 
boundary  layer  in  an  orthogonal  streamline  coordinate  system.  First-order  finite-difference 
methods  are  used  to  solve  the  resulting  boundary-layer  equations.  The  unknown  variables  are 
the  streannvisc  momentum  thickness,  the  shape  parameter,  and  the  streamline  slope  at  the  sur¬ 
face.  The  boundary-layer  calculation  method  is  combined  with  existing  geometrical  and  imiscld 
potential-flow  computer  codes  for  rotating  blades  to  form  an  efficient  turbulent  boundary-layer 
computer  code.  For  a  given  potential-flow  solution,  a  typical  boundary-layer  computation  re¬ 
quires  less  than  10  seconds  computer  time  on  the  Burroughs  7700  high-speed  computer. 
Boundary-layer  predictions  are  presented  for  several  rotating  blades.  Computed  results  are 
shown  to  be  in  agreement  wit!',  experimental  data  for  a  simple  rotating  body.  For  the  examples 
considered,  displacement  of  the  mid-chord  point  of  a  blade  from  a  straight  radial  line  is 
predicted  to  reduce  the  computed  values  of  local  skin  friction  coefficient,  with  an  estimated  in 
crease  in  overall  efficiency  of  about  one  percentage  point. 
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INTRODUCTION 


Turbulent  boundary  layers  occur  in  practice  on  the  rotating  blades  of  compressors,  turbines, 
helicopters,  or  propulsors.  C  entrifugal  and  Coriolis  forces  are  additional  terms  that  arise  in  the 
governing  equations  due  to  the  blade  rotation  relative  to  the  description  of  non-rotating  systems. 
These  forces,  along  with  the  pressure  and  viscous  forces,  result  in  a  complex  thrce-dimensioinn 
flow  configuration. 

Several  previous  studies  have  been  made  of  turbulent  boundary  layers  on  rotating  bodies. 
One  of  the  earliest  studies  was  by  von  Khrman***  in  1946  on  axisymmetrie  bodies.  Banks  and 
Ciadd<“\  in  (962,  investigated  both  laminar  and  turbulent  boundary  layers  on  rotating  blades.  \ 
review  of  much  of  the  recent  work  is  given  by  Meyne**)  (1972). 

In  1969.  C  hain  and  Head!**)  presented  experimental  and  theoretical  turbulent  boundary-laye* 
results  on  a  rotating  disk.  In  following  years,  they  extended  their  theory  to  include  the  flows 
over  a  rotating  circular  cylinder  n  an  axial  stream^)  and  a  rotating  axisvmmetric  bodv^T  \  1  i 
three  ot  Cham  and  Head's  investigations  used  the  momentum  integral*  "I  with  entrainment 
method  of  solution.  The  assumed  entrainment  function  and  the  exact  form  of  the  momentum  in¬ 
tegral  equations  varied  with  the  differing  geometries.  Ciood  agreement  was  l omul  between 
iheorebcal  pi  edict  ions  and  experimental  results  for  Cham  and  Head's  investigations. 

I  akshminai av ana,  el  ald‘9  presented  results  of  an  investigation  of  the  turbulent  bomulai  v 
lavet  on  a  rotating  helical  blade  in  19*72.  Again,  the  solution  was  based  on  the  momentum  in 
tegral  method.  I  he  governing  equations  were  simplified  by  the  assumption^  of  zero  pressme 
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gradient  and  the  use  of  the  orthogonal  streamline  coordinate  system.  Overall,  their  theoretical 
predictions  were  in  acceptable  agreement  with  their  experimental  results:  measured  values  of  I  he 
st reamwise  momentum  thickness  were  no  more  than  20  percent  less  than  predictions  at 
midradius  and  values  of  the  limiting  streamline  were  in  -+  2  degrees  agreement  with  predictions. 

Also  in  1972,  Meyne^)  presented  theoretical  and  experimental  results  lor  boundary-layer 
flow  on  rotating  wedged-shaped  surfaces.  Meyne  sought  a  generalized  similarity  solution  (reduc¬ 
ing  the  equations  to  an  ordinary  differential  equation  in  the  normal  direction)  and  considered 
both  laminar  and  urbulent  flow.  Detailed  results  are  shown  for  the  predicted  boundary-laser 
parameters  for  laminar  flow.  However,  his  equations  for  the  turbulent  flow  case  were  identical 
to  those  of  Banks  and  Ciadd^l  and  thus  no  new  information  was  obtained. 

The  present  work  extends  that  of  the  previous  investigators  by  providing  for  theoretical 
computations  of  the  turbulent  boundary  layer  that  develops  about  more  complex  geometries.  In 
particular,  the  current  method  was  designed  to  compute  the  turbulent  boundary  layer  on  the  sur¬ 
faces  of  rotating  blades  with  thickness  and  loading  in  an  axial  stream  operating  at  0. 7-radius 
Reynolds  numbers  greater  than  about  4  \  10^  for  which  fully  turbulent  flow  is  expected. 
Although  not  considered  explicitly,  the  present  calculation  method  could  be  modified  slightly  to 
perform  the  laminar  calculation  as  well. 

The  three-dimensional  turbulent  boundary  layer  in  streamline  coordinates  (4.  >].  4)  is  govern 
ed  by  three  equations.  These  are  the  4  and  q  momentum  integral  equations  and  the  continuity 
equation.  These  equations  are  the  same  as  those  given  by  SmitMhh  -md  by  Oimpstv  and 
Head! 1 1  *  with  one  additional  term  (to  account  for  the  rotation)  in  both  of  the  momentum  in¬ 
tegral  equations.  In  the  integral  method,  the  equations  of  motion  are  integrated  across  ihe  boun¬ 
dary  layer,  producing  first-order  partial  differential  equations.  Tmpirical  information  is  incot- 
poruted  in  the  equations  to  account  for  Reynolds  stresses.  In  the  numerical  method,  diiect  use  to¬ 
rnado  of  the  streamwise  momentum  thickness  fljj  and  the  shape  parameter  H,  the  two 
parameters  which  most  readily  characterize  the  behavior  of  turbulent  bouiularv  lavcisH2)_  | 


use  of  streamline  coordinates  simplifies  both  the  boundary  conditions  and  the  governing  equa¬ 
tions  but  requires  the  definition  of  a  new  coordinate  system  based  on  the  How  field  and  blade 
geometry. 

Several  assumptions,  which  are  discussed  in  a  later  section  of  this  paper,  are  made  concern¬ 
ing  the  behavior  of  the  How.  These  assumptions  reduce  the  number  of  unknown  quantities  in 
the  three  governing  equations  to  three  variables.  These  are  the  streamwise  momentum  thickness 
0 1  j,  the  shape  parameter  H,  and  the  tangent  of  ft,  the  angle  between  the  direction  of  the  wall 
friction  vector  (or  wall  streamline)  and  the  inviscid  streamline. 

When  computing  the  turbulent  boundary  layer  on  a  rotating  blade  with  thickness  and  lift,  a 
potential-flow  computer  code  1 1 developed  at  DTNSRDC  by  Brocken  is  used.  This  code  com¬ 
putes  the  potential  flow  velocities,  slopes  of  the  streamlines  and  crosswise-normal  lines,  metric 
coefficients,  and  the  angular  velocity  component  normal  to  the  surface.  All  values  are  given  in 
non-dimensional  coordinates  (xc,  \K)  where  xL.  is  the  chord  fraction,  measured  from  the  leading 
edge  and  xR  is  the  fraction  of  radius,  measured  from  the  axis  of  rotation. 

Streamlines  and  crosswise-normal  lines  are  determined  from  t he  given  slopes  using  an  Lulor- 
predictor/trape/oidal-corrector  method*^.  Potential-How  velocities,  metric  coefficients,  etc.  arc 
interpolated  in  the  streamline  coordinate  system  using  a  spline-under-tension  method!  I-Sl.  A  first- 
order  finite-difference  approximation!  is  used  to  solve  the  system  for  the  three  basic  unknown 
quantities. 

No  interaction  with  the  hub  is  considered  in  the  present  analysis  and  the  nature  of  the  flow 
near  the  tip  is  assumed  to  be  adequately  predicted  by  the  potential- flow  model. 

Results  of  the  turbulent  boundary-layer  calculation  are  presented  for  several  rotating 
geometries.  The  first  body  studied  is  the  rotating,  weakly  helical,  segment  investigated  ex¬ 
perimentally  and  analytically  by  I  akshminarayana,  el  al.^h  Computed  results  are  in  agreement 
with  the  experimental  data.  The  remaining  geometries  are  described  unalylicnIK  with  no  previous 
experimental  or  analytical  boundary-layer  data  available.  I  lie  effect  of  blade  warp  and  skew  on 
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the  predicted  values  of  the  local  skin-friction  coefficient  is  investigated  qualitatively  using  these 
surfaces.  Blade  warp  is  defined  as  the  angular  displacement  of  the  midchord  point  from  the 
reference  line  in  the  plane  of  rotation  and  blade  skew'  is  defined  as  the  displacement  of  the  mid¬ 
chord  point  along  the  pitch  helix  at  constant  radius. 


COORDINATE  REFERENCE  FRAMES 

The  geometry  and  boundary-layer  characteristics  of  rotating  blades  ate  described  by  two 
coordinate  systems  fixed  in  the  rotating  surface.  Brocket!*1-)  gives  the  geometry  and  potential 
flow  characteristics  in  terms  of  the  non-dimensional  coordinates  xL.  and  \^.  see  figure  1.  I  hese 
are  uon-orthogonal,  helical  coordinates  where  \c  is  the  fraction  of  chord  as  measured  from  the 
leading  edge  and  x^  is  the  fraction  of  radius  as  measured  from  the  axis  of  rotation. 


figure  I.  Blade  Surface  Cieometix 


The  turbulent  boundary-layer  scheme  uses  an  orthogonal,  curvilinear  coordinate  system  bas¬ 
ed  on  the  potential-flow  streamlines,  but  with  explicit  dependence  on  the  xc  and  xR  coordinates. 
The  use  of  the  orthogonal  coordinates  simplifies  the  boundary  conditions  and,  hence,  the 
boundary-layer  equations  governing  the  turbulent  now.  The  orthogonal  streamline  coordinates 
on  the  surface  are  4,  which  varies  along  streamlines,  and  q,  which  is  normal  to  the  streamline. 
The  4  coordinate  measures  distance  along  the  surface  normal. 

In  the  orthogonal,  curvilinear  coordinates,  the  distance  between  lines  of  constant  4  or  rj 
varies  from  point  to  point  on  the  surface.  Metric  coefficients,  denoted  by  h|,  I12,  and  hj  for  the 
4.  q,  and  4  coordinates,  respectively,  are  required  to  1  arrclaic  surface  distances  to  increments  of 
the  streamline  coordinates.  Since  the  boundary-layer  thickness  is  assumed  to  be  small  and  the  4 
coordinate  is  assumed  to  measure  distance  along  the  surface  normal,  the  metric  can  be  set 
equal  to  unity  with  no  loss  of  generality.  A  general  element  of  length  ds  in  streamline  coor¬ 
dinates  is  then  given  by 

ds-  =-  h,  -dO  4  h?  -dr]-  4  d£-  (I) 

where  h|  and  hi  are  functions  of  4  and  r]  only. 

BASIC  EQUATIONS 
Momentum  liquations 

The  continuity  and  momentum  equations  for  laminar  flow  in  a  totaling  system  are  readily 
available  in  the  literature,  e.g.,  Moore* *  'h  The  turbulent  boundary-layer  equations  are  the 
laminar  equations  plus  additional  terms  due  to  the  Reynolds  stresses.  In  streamline  coordinates, 
the  equations  for  turbulent  flow  are 
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where  u,  v,  and  w  are  the  time-mean  velocity  components  in  the  4.  >1  and  4  directions,  respectively, 
K,  is  the  geodesic  curvature  of  lines  of  constant  |, 


K]  =  — L  — 

h|h2 


k->  is  the  geodesic  curvature  of  lines  of  constant  r), 
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U  is  the  speed  along  the  streamline  at  the  outer  edge  of  the  boundary  layer, 

[a  is  the  viscosity  of  the  fluid, 
q  is  the  density  of  the  fluid, 

oji  is  the  angular  velocity  component  in  the  4-direction  (n -direct ion ) 

go}  =  co  ■  n  =  -  cu(j_  '  n ) 

u'w'  and  v'w'  are  the  time  averaged  products  of  the  turbulent  fluctuating  velocity  components 
(i.e.,  components  of  the  Reynolds  stresses). 

The  integral  boundary-layer  equations  are  obtained  by  integrating  the  momentum  equations 
with  respect  to  4  from  the  surface  of  the  blade  (4  =  0)  to  the  edge  of  the  boundary-layer  (4—3). 
The  continuity  equation  is  used  to  eliminate  all  terms  involving  the  velocity  component  \\,  nor¬ 
mal  to  the  blade  surface.  The  resulting  equations^)  are: 


4 -moment uni  integral: 


1  00ii  13012 

—  -  + - —  +  ho(0 1 1  ~  #22)  +  K  |  (202 1  -  d2) 

hj  34  h2  3r] 

2021-92  dU  20]  |+9,  9U  2w3  ,  \_r. 

Uh2  9  rj  +  Uh,  94  U  2  2  f| 


(3a) 
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rj-momentum  integral:  —  ^21  +  —  ^22  +  20,1  (  —  —  +  K? 

h,  di  h2  dr,  21  V  Uhj  04  2 

2uj  2022  3U  ] 

+  —  <’1  +TSIT,  -  ^<01,-922^,1  =Tcr: 


13b) 


In  equations  3a  and  3b,  0jj,  012,  02i,  022’  and  d2  are  the  momentum  and  displacement 
thicknesses  defined,  respectively,  by  succeeding  equations  5a  through  51  and 

d  v 


C>l  = 


du  — — 

“at  "euw 

yeu2 


and  C}7  = 


w  —  -  e'  w 
34 


yeu2 


(4) 


arc  the  turbulent  wall  skm  friction  coefficient  components  in  the  4  and  r,  directions,  respectively. 
I  he  momentum  and  displacement  thicknesses  are  defined  by: 


8"  - !('  -b)bdt 


(5a) 


'“-'('-irfe* 


(5b) 


„  r  UV 

021  =  -;-T  d4 
0 

d2  ♦  0 1 


(5c) 


d  v2 


e22  =  -/  ry 

O  U“ 


( 5d ) 


--jjjdi 


(5c) 


d2  =  -  /—  d£ 
J0U  4 


(50 


The  boundary  conditions  are  u  =  v  =  0  at  4  -  0  and  u  U,  \  0  at  4  =  0,  the  edge  ol 

the  boundary-layer. 
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The  two  momentum  integral  equations  are  insufficient  to  define  the  boundary-layer 
development.  An  additional  equation  (the  entrainment  equation)  will  be  introduced  in  the  next 
section  and  several  assumptions  (for  C^,  Cf2,  and  mean  velocity  profile  shapes)  applied  in  order 
to  reduce  the  number  of  unknown  quantities  to  three.  These  three  quantities  are  the  streamwise 
momentum  thickness,  0n;  the  shape  factor  parameter,  H;  and  the  tangent  of  the  angle  between 
surface  streamline  and  external  streamline  directions,  tan/3. 

Entrainment  Equation 

The  three-dimensional  entrainment  equation  of  MyringW,  which  has  also  been  used  by 
Cumptsy  and  Head^1^  and  von  Kerczek  and  Langan*18),  is  the  additional  equation  used  here  to 
describe  the  turbulent  boundary-layer  How  over  the  blade  surface.  This  equation  describes  the 
growth  of  the  boundary  layer  with  time  as  the  turbulent  mixing  process  spreads  the  general  mo¬ 
tion  of  turbulent  flow!19).  The  continuity  equation  is  integrated  with  respect  to  4  through  the 
boundary  layer  to  obtain  (in  streamline  coordinates)  the  entrainment  equation: 


1  d(d-dj)  1  d&2  +  |  1 


hj  34  h2  dn 
1  3E 


Uhl  n  +  «-'■* 


^  ^  I^U  dr\  + 


,  1  3d 

K')  =  V,  JT 


w(c5) 


(6a) 


h|  di  U 

The  right-hand  side  of  equation  (6a)  is  a  measure  of  the  rate  at  which  fluid  enters  the  boundary 
layer.  Denoting  this  quantity  by  F,  where  F  is  a  dimensionless  entrainment  function,  the  entrain¬ 
ment  equation  becomes: 


1  d(d-d])  1  dd2 


1  3tH  !  1  3U  1/  \  J  r  r  V  .  /  1  ^  _  IV 

(  ui^  aT  *  k->  - *  (  HJU  XT  *  K|>  '  ' 


1  3U 


H 


(6h) 


Details  of  the  entrainment  function  are  given  in  a  later  see 


■linn  nt’  l  hiv.  nnrw'r 


ADDITIONAL  RELATIONSHIPS 

Further  assumptions  are  needed  in  order  to  solve  the  two  momentum  integral  equations  (3a 
and  3b)  and  the  entrainment  equation  (6a). 
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Skin  Friction 


The  turbulent  skin  friction  terms  Cf ^  and  Cf2  on  the  right-hand  side  of  the  £  and  y]  momen¬ 
tum  integral  equations,  respectively,  are  obtained  from  empirical  data  for  the  skin  friction  coef¬ 
ficient  Cf  in  the  local  inviscid  streamline  direction^8)  by 


Cf,  =  Cf 

(7a) 

and 

Cf2  =  Cf  tan  ft 

(7b) 

where 

3  v 

v  34 

tan/3  =  lmi  —  =  lim  du 

4-o  4—o 

(7c) 

By  assuming  that  the  component  of  boundary-layer  flow  in  the  inviscid  streamline  direction 
is  approximately  two-dimensional,  the  skin  friction  coefficient  Cf  can  be  specified  using  a  two¬ 


dimensional  skin  friction  formula.  The  skin  friction  relationship  given  by  Thompson^))  ancj  ap¬ 
proximated  by  Head  and  Patel^D  is  used  in  the  present  calculation  method.  The  skin  friction 


coefficient  is  approximated  by  the  analytic  relation 

Cf  -  exp  (aH  +  b)  (8a) 

where  a  -  0.019521  -  0.386768c  +  0,028345c2  -  0.000701c3  (8b) 

b  -  0.191511  -  0.834891c  4-  0.062588c2  -  0.001953c3  (8c) 

c  -  ,nR011  (Sd) 


U0|| 

Re,,  — 

U  =  speed  at  the  edge  of  boundary  layer  (8g) 

v  =  kinematic  viscosity  of  water  ( Sh ) 

6 1 1  -  the  dimensional  streamline  component  of  the  momentum  thickness.  ( Si ) 

Both  coefficients  a  and  b  decrease  as  #j|  increases.  Thus,  for  decreases  or  “small"  increases  in 
the  shape  factor  H,  the  skin  friction  coefficient  C(  decreases  as  #jj  increases. 
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Velocity  Profiles 


Assumptions  involving  the  formulation  of  the  streamwise  and  crosswise  velocity  profiles  are 
made  in  order  to  relate  the  crosswise  momentum  and  displacement  thicknesses  d\2<  02l>  e22>  and 
d2  to  the  streamwise  momentum  and  displacement  thicknesses  0ji  and  A j.  Extreme  accuracy  is 
not  necessary  in  the  representation  of  these  velocity  profiles  since  they  are  used  only  to  relate  in¬ 
tegrated  values  of  the  crosswise  and  streamwise  momentum  and  displacement  thicknesses. 

The  approximation  tor  the  streamwise  velocity  profile  is  given  by  Cooke*2-)  and  by 
Patel*-^.  This  formulation,  which  approximates  the  experimental  swept-wing  boundary-layer 
data  of  W  allace*24),  is 


H-l 


(9a) 


The  crosswise  velocity  profiles  are  assumed  to  be  related  to  the  streamwise  velocity  profile^ 


by  empirical  formulation  due  to  Mager*2^.  This  approximation,  also  assumed  by  Cooke*22*, 
C  ham  and  Head*4*,  and  von  kerc/ek  and  Eangati*^),  is  given  by  the  relation 


U 


TJ  * 1  •  7  )‘ tan/? 


(9b) 


This  approximation  is  assumed  to  be  adequate  as  long  as  it  is  not  applied  to  crossflow  reversal. 

An  approximation  for  the  boundary-layer  thickness  6  is  found  by  substituting  the  assumed 
form  of  the  streamwise  velocity  profile,  equation  (9a),  into  the  definition  of  the  streamwise 
momentum  thickness  0j|,  equation  (5a).  Integration  and  simplification  give  the  approximation 


'>-<’11 

C  rosswise  I  hicknesses 

The  expressions  for  streamwise  velocity  profile,  equation  (9a),  and  the  crosswise  velocitx 
profile,  equation  (9b)  are  substituted  into  the  definitions  of  the  streamwise  displacement 
thickness  and  the  crosswise  momentum  and  displacement  thicknesses,  equations  (5b)  through 


1 


(5f),  and  integrated  across  the  boundary  layer.  After  considerable  manipulations,  the  following 

expressions  result: 

d\2 

=  0U  MI2 

( 10a) 

021 

=  6n  m2, 

(10b) 

e22 

=  011  m22 

(10c) 

*1 

=  011  Dj 

(10d) 

<*2 

=  011  d2 

(10c) 

4 

1 

4  2  v 

where  M ]2  =-(G 

+  H)  tan/?  ( 

H  +  l"+  H  +  2  1  H  +  3  H  +  5  ' 

(11a) 

m2,  = 

=  -(G  +  H)  tan/? 

(rr 

2 

H  + 1  + 

1  ) 

H  +  2  > 

(lib) 

M22  = 

=  -(G  +  H)  tan2/? 

(1 

4 

H  +  1  + 

6  .  4  +  1  ) 

H  +  2  H  +  3  H+4  ' 

(11c) 

D,  = 

H 

(lid) 

D2  = 

-(G  +  H)  tan/?  ( 

2 

H  +  1 

4 

H  +  3  + 

2  ) 

H  +  5  ! 

(lie) 

Entrainment  Function 

The  entrainment  function  F  is  assumed  to  be  determined  uniquely  by  the  circumferential 
velocity  profile  in  the  same  way  as  proposed  by  Head*19*  for  a  two-dimensional  boundary  layer. 
However,  C  ham  and  Head<4)  found  that  if  the  entrainment  function  is  reduced  to  approximately 
two-thirds  of  the  Hat  plate  value,  then  agreement  between  calculated  turbulent  boundary-layer 
development  on  a  rotating  disk  and  experiments  was  excellent.  Incorporating  the  results  of 
Cham  and  Head  into  the  entrainment  function  presented  by  Cireen,  el  al.(-^  gives  the  following 
universal  expression  for  the  function  F: 

F  =  0.025  (TL-1) 
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G 


2H 

H-l 


(12b) 


H  =  the  shape  factor  parameter,  H  =  -J-  (12c) 

The  relationship  between  the  boundary-layer  thickness  6  and  the  streamwise  moment um 
thickness  0]j  given  in  equation  (9c)  can  be  rewritten  using  the  definition  of  G  as 

6  -  0j|  (G  +  H)  ( 1 2d) 


The  three  basic  equations  which  characterize  the  turbulent  How  over  the  surface  of  the 
rotating  segment  (the  momentum  integral  equations  (3a)  and  (3b)  and  the  entrainment  equation 
(6b))  are  now  expressed  in  terms  of  the  three  basic  unknown  quantities:  the  streamwise  momen¬ 
tum  thickness  6 jj,  the  shape  factor  parameter  H,  and  the  tangent  of  the  angle  between  surlace 
streamline  and  external  streamline  directions,  t  -  tan/3. 


NUY1TR1C  AI  ANA1  V  SIS 

In  order  to  obtain  a  more  compact  form  of  the  three  basic  equations,  (3a),  (3b),  and  (6b), 
the  expressions  for  the  crosswise  momentum  and  displacement  thicknesses  and  the  streamwise 
displacement  thickness,  equations  (l()a)  through  (l()e),  are  substituted  into  each  of  these  equa¬ 
tions.  After  some  manipulation,  the  single  hyperbolic  equation  in  vector  form  becomes 


a  w  a  w 

A(W  (In))  —  +  B(W(i,r)))  —  =  C(W) 

dk  or] 

[  he  unknown  vector  W  is  given  by 
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and  the  known  coefficient  matrices  are 


i  d  Mi  1 

T7  ~dT 


9 1 1  d  M  2 ) 

ITT  "~3TT 


Oj  1  a  (G 

+  H-D|) 

0i  1 

hi 

at 

hl 

a,,  a m  12 

01! 

a  m  ,  2 

h2  at 

h2 

3H 

9  j  j  a  m  n 

el  1 

a  M22 

b2  at 

b2 

3H 

a,,  a d2 

0ii 

3[)^ 

‘  h2  at" 

h  2 

Tff 

c'i  =■-  tc'i,  -  K:0i  1  O-M::)  -  Mu  <2M21-d2>  •  9,1 


0ii(2*i)|)  au 


U  h,  34  U  0,1  D: 

( 2  -  3C  t2  ■  ~0iim2i  (  "ay  4  K2 )  ■  _[^eiiI)i 
+  K  j  $i  |(  1  Nb 2  *  ) 

f  2  1  -  01 1  <(i  •  H  - *  k2  )  •  0||l->2  ( 


20| |\i22  a i 
h  2  l  d  r\ 


1  ai_ 

h  ■»  l  d  n 
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Since  equation  (13)  is  hyperbolic,  the  entire  solution  between  asymptotes  could  be  determin¬ 


ed  once  initial  conditions  were  specified  along  these  limiting  lines.  Initial  conditions,  however, 
are  not  readily  available  for  the  geometries  of  interest.  Therefore,  a  solution  procedure  is 
adopted  which  yields  results  only  along  one  side  of  the  bisector  of  the  asymptotes  and  requires 
initial  conditions  along  only  one  asymptote. 

A  first-order  finite  difference  approximation  is  used  to  solve  equation  (13)  in  the  streamline 
coordinates  (4,  >7).  The  grid  spacing  is  given  by  A4  in  the  st teamwise  direction  and  Arj  in  the 
crosswise  direction.  Arclengths  are  approximated  by  the  chordlength;  hence,  reasonably  small  in¬ 
crements  are  required.  Superscript  i  refers  to  the  crosswise  (p)  direction  and  superscript  i  refers 
to  the  stream  wise  (£)  direction,  forward  differencing  in  4  yields 

( a  w\  i  w-i  *  1  -  vvj 

V  j  ~  hi  ( !5a) 

and  backward  differencing  in  n  yields 

/ d\\  \  1  \V‘  -  VV'-’ 

\drl  )  ^  ( I5h| 

II  the  \  allies  ot'  H,  and  tan  ft  are  known  at  the  grid  node  (i,j),  the  values  of  these  three 
variables  can  be  found  at  the  node  (ij  4  1)  by  (he  first  order  finite  difference  approximation 


AM 


j'VVM 


*  1  -  vvi.A 
hi  } 


BU 


(W'ij  ■  VV'  l-A 

\~An~J 


(i.i 


It  the  A1*1  matrix  is  nonsingular,  equation  (16)  can  be  solved  for  \\  '-i 


(  16) 


VV'-j  ' 

1  - 

£  [(AM) 
An 

1  B'-i  ( w  '•> 

-  W'-'-jfl-t  W>-j 

<r» 

where 

hi 

1  -  V 

( 1  Sal 

and 

An 

£»i'  -  n1'1 

<1M'. 
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The  streamline  grid  is  shown  on  a  typical  sector  rotating  at  a  constant  angular  speed  ca  in 
Figure  2.  As  experiments  indicated,  streamlines  (lines  along  which  £  varies)  are  shown  to  ap¬ 
proximate  constant  radius  lines  and  the  crossflow  lines  (lines  along  which  rj  '  arics)  are  shown  to 
approximate  radial  lines.  Streamlines  are  indicated  by  the  variable  j  and  crossflow  lines  by  the 
variable  i.  The  arelengths  A£  and  A  rj  are  shown  to  vary  from  point  to  point  on  the  sector  sur¬ 
face.  Grid  nodes,  which  are  defined  as  intersections  of  streamlines  and  crossflow  lines,  are  label¬ 
ed  i,i  for  i  -  1  to  N,  j  -  1  to  M  where  N  -  number  of  streamlines  and  M  number  of 
crossflow  lines. 


CROSSFLOW  LINES 

1  igurc  2.  Streamline  Grid  on  a  1  \ picul  kotating  Sector 


1  igurc  3  shows  the  calculation  sequence  for  the  (urlmlcni  boundan  laver  on  a  lvpical  see 
.oi  1  or  convenience,  the  grid  nodes  used  in  the  calculation  procedure  arc  numbered  sequent  in  !l\ 
from  the  number  1.  It  is  assumed  that  starting  values  of  6 jj,  11,  and  t  tan/5  arc  given  along  a 
crosstlow  line,  points  numbered  1-7.  Turbulent  boundary-layer  characteristics  are  solved  along 
successive  crossflow  lines  as  indicated  by  the  increasing  numbers  in  Figure  3.  I  or  example,  dat  i 
at  points  j  and  2  are  used  to  obtain  predictions  at  point  K,  data  at  points  2  and  3  aic  used  foi 
predictions  at  point  9,  etc.  As  can  be  seen  from  this  simplified  example  of  the  calculation  pio 
ccdure .  a  line  grid  is  necessary  near  the  hub  of  the  blade  in  order  to  obtain  results  near  the  hail 
mg  edge. 


A  blade  is  typically  defined  by  a  grid  containing  about  250  nodes.  Experience  to  date  in¬ 
dicates  that  grid  nodes  should  be  distributed  on  the  blade  with  a  greater  concentration  near  the 
inner  radii  and  leading  edge;  and  arclengths  A£  and  Arj  should  increase  uniformly  in  the  direc¬ 
tions  of  the  trailing  edge  and  tip  of  the  blade.  Using  the  Burroughs  7700  high-speed  computer, 
the  boundary -layer  computation  for  this  grid  requires  less  than  10  seconds  computer  Pme  at  a 
cost  of  less  than  one  dollar. 


figure  3.  Points  at  Which  Calculations  Performed 

COMPUTATIONAL  RESULTS  AND  DISCUSSION 

The  computer  code  for  calculating  the  three-dimensional  turbulent  boundary  layer  has  been 
evaluated  using  several  geometries  for  which  analytical  and  experimental  data  are  available. 
Among  these  are  a  flat  plate,  a  two-dimensional  foil,  a  rotating  disk,  and  a  rotating  helical 
blade.  Excellent  agreement  with  test  data  was  tound  in  each  case.  Lite  turbulent  boundary  lav  or 
was  computed  for  the  rotating  helical  blade  assuming  both  zero  and  nonzero  crossflow.  A‘  was 
anticipated,  results  of  these  two  computations  showed  agreement  within  15  percent.  Results  arc 
presented  for  only  the  more  complex  of  these  geometries,  the  rotating  helical  blade  with  non/ct 
crossflow.  The  turbulent  boundary-layer  development  over  the  surface  of  Blade  1,  a  sector  of 
300  degrees  included  angle,  was  investigated  experimentally  and  analvticallv  bv 
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Lakshminarayana,  et  al.W.  The  small  pitch  of  the  sector  was  not  included  in  their  analysis,  nor 
is  it  included  in  the  present  investigation. 

The  effect  of  blade  warp  and  skew  is  investigated  using  three  blades:  an  unskewed  Blade 
2A,  a  72  degree  warped  Blade  2B  and  a  72  degree  skewed  Blade  2C.  Blades  2  are  analytically- 
defined  lifting  surfaces  for  which  no  previous  experimental  or  analytical  turbulent  boundary- 
layer  data  exist.  Details  of  the  geometry  and  potential  flow  calculations  are  given  in  Reference 
13.  Initial  conditions  for  this  geometry  are  assumed  to  be  flat  plate  and  rotating  disk  values  for 
turbulent  flow.  The  computed  local  skin  friction  coefficients  are  compared  for  these  blades.  At 
present,  only  the  skin  friction  coefficients  are  computed.  In  the  future,  it  is  anticipated  that  the 
force  and  moment  values  will  also  be  computed,  allowing  a  calculation  of  the  effect  of  viscous 
drag  on  overall  performance. 

Blade  1 

The  turbulent  boundary-layer  characteristics  were  investigated  initially  on  a  single  non¬ 
lifting  rotating  helical  sector  of  92.6  cm  (36.6  in.)  diameter  with  a  300  degree  included  angle, 
hence  of  a  large  chord  length.  Both  experimental  and  analytical  data  are  available  from 
Lakshminarayana,  et  ald^  for  comparison. 

1  he  pitch  ot  the  helical  surface  studied  by  l  akshminarayana,  et  ald^l  was  25.4  cm  (10  in¬ 
ches),  tor  a  pitch  diameter  ratio  ot  0.273.  Muis,  the  geometry  can  be  approximated  by  a  flat  cir¬ 
cular  disk  with  a  leading  and  trailing  edge  300  degrees  apart.  The  blade  was  mounted  on  a  22.86 
cm  (0.75  ft)  radius  hub  and  the  maximum  radius  of  the  blade  was  46.48  cm  (1.525  ft).  The  blade 
was  enclosed  in  a  housing  wall  with  a  0.20  cm  (0.08  in.)  clearance  between  the  blade  tip  and  the 
wall.  The  frec-siream  onset  flow  was  zero  and  the  rotational  speed  of  the  blade,  denoted  to,  re¬ 
mained  constant  at  47  rad/s  (450  rev  min).  The  geometry  of  the  sector  is  shown  in  f  igure  4. 

The  theoretical  analysis  of  Lakshminarayana,  et  al.<^  is  based  on  a  rotating  system  describ¬ 
ed  with  cylindrical  coordinates  (r,  6 ,  /)  where  r  varies  in  the  radial  direction.  8  varies  in  the 
chordwise  direction  (i.e.,  circumferential  direction)  and  /  varies  along  the  surface  normal  (axial 
direction).  These  coordinates  are  indicated  in  f  igure  4  and  are  used  when  comparing  the  ex¬ 
perimental  data  with  the  results  of  the  present  theory. 
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f  igure  4.  Geometry  and  Coordinates  of  Rotating  Helical  Sector 
The  potential-flow  surface  velocity  for  Blade  1  is  simply  cur  along  the  chordlines.  The 
streamlines  arc  lines  of  constant  radius  along  which  neither  the  potential  flow  velocity  nor  the 
metric  coefficient  h^  R  varies.  The  momentum  integral  and  entrainment  equations  for  this 
sector  are  thus  simplified  to: 

$  momentum  integral : 


r)4  h: 
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l  lie  metric  coefficients  h  j  and  h2  are  the  chord  length  and  maximum  radius,  respectively  and  4 
and  y]  are  the  changes  in  the  fraction  of  chord  and  fraction  of  radius,  respectively.  Substitution 


of  these  values  into  equations  (19a)  and  (19b)  yields  the  momentum  integral  equations  given  by 
Lakshminarayana,  et  al.W.  In  their  analytical  study,  only  the  streamwise  momentum  thickness 
$1  [  and  the  limiting  streamline  angle  ft  varied.  A  general  form  of  the  solution  was  assumed  (0|  j 
was  assumed  to  vary  with  radius  and  chord  and  tan/?  was  assumed  to  vary  with  0 ]  j  only)  and  the 
governing  momentum  integral  equations  solved  numerically  using  the  fourth-order  Runge-kutta 
method.  In  reference  (9),  it  was  assumed  that  the  streamwise  velocity  profile  varies  as  the  1  7 
power  of  (c  d)  and  a  simple  friction  formula  was  selected.  The  present  method  also  uses  the  en- 
iiainment  equation  and  a  different  solution  method  based  on  assumed  entrainment  and  skin  fric¬ 
tion  relationships  and  assumed  velocity  profiles.  Thus,  some  differences  occur  in  the  predicted 
flow  parameters. 

!  akshminarayana,  et  al.M  presented  experimental  and  theoretical  results  at  a  rotational 
speed  of  4"  rad  s  (450  res  min).  The  kinematic  viscosity  v  equaled  1.486  x  l(Ms  m~/\  (160  \  l(M' 
ft*-  sec)  which  corresponded  to  air  at  20  C  (68rT).  These  conditions,  for  a  Reynolds  number 
based  on  tip  radius  and  tip  rotational  speed  of  7  \  10s  (or  a  0.7-radius  Reynolds  number  of  2.5 
\  10^  (see  next  section  for  definitions)),  are  used  in  the  evaluation  of  the  present  method. 

In  the  numerical  solution  procedure,  I. akshminarayana,  et  ald^  performed  calculations  ot 
the  laminar  boundary  layer  until  the  flat-plate  critical  Reynolds  number,  Rv  u)\~6  v  5  \ 

!0-\  was  reached.  The  salues  of  the  streamwise  momentum  thickness  f?jj  and  limiting  streamline 
angle  6  weie  then  used  as  initial  conditions  for  the  turbulent  calculations,  blow  \ isiiali/ation  ex¬ 
periments  using  a  sublimation  technique  supported  the  use  of  the  flat  plate  critical  Res  Molds 
number  to  determine  the  transition  point. 

With  the  present  scheme,  calculations  began  at  the  location  6  OM33  radians  and  xK 
i  R  0.492,  the  grid  node  location  closest  to  the  experimentally-determined  transition  location. 
Initial  conditions  for  the  streamwise  momentum  thickness  0jj  and  the  tangent  of  the  limning 
streamline  /?  were  obtained  from  the  experimental  data  of  I  .akshminarayana,  et  ald^h  but  no 
radial  variations  were  considered.  Initial  values  were  taken  to  be  6 j|  0.24  mm  (0.0008  ft), 
tan/I  0.2,  and  H  1.5  at  each  radial  location  along  the  chordsvise  position  0  0.733 

i  adians. 
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Figures  5  through  7  present  the  results  of  the  turbulent  boundary-layer  calculation  over  the 
surface  of  Blade  1.  Both  the  experimental  and  analytical  results  of  Lakshminarayana,  et  ai.^1 
and  Oat  plate  predictions  (adjusted  for  the  region  of  laminar  How  prior  to  0  =  0.733  radians) 
are  indicated  in  the  figures  for  comparison.  The  variation  in  the  stream  wise  momentum 
thickness  is  given  in  Figures  5a  through  5c  at  the  radial  locations  xR  ~  r/R  -  0.72,  0.82, 
0.93,  respectively.  The  momentum  thickness  is  predicted  to  be  less  than  the  flat-plate  value  given 
by  Prandll  approximation 1  (0 j  j  0.036.x  Rv  1  with  U  -  cur  and  \  rfl),  and  less 
than  the  prediction  of  Lakshminarayana,  et  al.^K  Use  of  the  present  theory  gives  satisfactory 
agreement  with  the  experimental  results. 

The  computed  variation  of  the  limiting  streamline  angle  ft  is  shown  in  f  igure  6a  for  r  R 
0.74  and  in  Figure  6b  for  r  R  0.965.  The  peak  in  the  surface  flow  angle  at  the  leading  edge  is 
the  result  of  laminar  flow.  The  present  calculations  start  at  6  -  0.733  radians  (42  degrees),  the 
measured  point  of  transition  at  i  R  0.72.  These  figures  show  that  the  currently-applied 
analytical  model  predicts  greater  variations  in  the  limiting  streamline  angle  than  observed  for  the 
measured  data.  The  laminar  sublayer  is  not  included  in  the  present  calculations  and  mav  be  im¬ 
portant  i cm  mis  variable.  Values  near  the  tip  may  be  influenced  by  the  annular  wall. 

The  distribution  of  the  local  skin  t  net  ion  coefficient  for  varying  radial  locations  is  given  in 
figure  Ta  at  (9  2.64  radians  and  figure  7b  at  0  4.76  radians.  Flat  plate  predictions  were 

made  using  the  SclioenherH-^l  curve  for  local  skin  friction,  (  t  0.310  |ln~(2Rfp  f  21n(2R^)J. 
figure  ~  shows  that  (lie  measured  values  of  the  local  skin  friction  coefficient  are  huger  than  die 
flat  plate  predictions  and  smaller  than  the  experimental  predictions. 

Blades  2 

In  this  section,  icsiilts  are  presented  for  the  computed  boundary-layer  charaetei isties  on 
three  artalv ( icallv  defined  blades*®  'L  an  imuarped  propuhor,  Blade  2A  (similar  to  NNRIX 
Model  4 4 S l ) ,  a  lughlv  watped  piopukor.  Blade  2B  (similar  to  \SRIX  Model  441>S),  and  a 

I  Prandll *s  appioximation  was  adjusted  for  the  region  of  laminar  flow  h>  matching  the  values  ot 
laminar  and  turbulent  momentum  thickness  #|j  at  the  location  0  0.733  radians. 
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Figure  5.  Variation  of  Momentum  Thickness  0jj  on  Rotating  Helical  BlaJc. 
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I  itMiie  ”,  Variation  of  Skin  Iriction  Coefficient  Cf  on  Kotatinu  Helical  If  lade. 
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The  geometry  of  the  warped  blade  is  given  in  Table  1.  The  geometry  of  the  unwarped  blade 


is  identical  with  the  exception  that  the  skew  angles  are  zero.  The  blades  rotate  about  a  cylin¬ 


drical  hub  at  a  constant  angular  velocity  cu  in  a  constant  axial  onset  stream  of  speed  V. 


TABLE  1  —  GEOMETRY  OF  ROTATING  BLADE  2B 


Fraction 
of  Radius 

Chord 

Pitch 

Skew  Angle 
(radians) 

Thickness 

Diameter 

Diameter 

Chord 

0.20000 

0.16500 

1.49500 

0.00000 

0.24000 

0.25000 

0.19700 

1.53500 

0.07854 

0.19800 

0.30000 

0.22900 

1.54400 

0.15708 

0.15610 

0.40000 

0.27500 

1.48700 

0.31416 

0.10680 

0.50000 

0.31200 

1.38900 

0.47124 

0.07680 

0.60000 

0.33700 

1.27700 

0.62832 

0.05650 

0.70000 

0.34700 

1.15600 

0.78540 

0.04210 

0.80000 

0.33400 

1 .04200 

0.94248 

0.03140 

0.90000 

0.28000 

0.93400 

1 .09956 

0.02460 

0.95000 

0.24000 

0.87700 

1.17810 

0.02330 

1.00000 

0.00000 

0.81200 

1.25664 

0.02460 

Fraction 

Thickness 

Chord  wise 

of  Chord 

Offset 

Loading 

0.000000 

0.000000 

0.000000 

0.007596 

0.086824 

0.912003 

0.030154 

0.171010 

0.992645 

0.066987 

0.250000 

1.040899 

0.116978 

0,321394 

1.074103 

0.178606 

0.383022 

1.097916 

0.250000 

0.433013 

1.114882 

0.328990 

0.469846 

1.126317 

0.413176 

0.492404 

1.132938 

0.500000 

0.500000 

1.135109 

0.586824 

0.492202 

1.135109 

0.671010 

0.467219 

1.135109 

0.750000 

0.421900 

1.103578 

0.821349 

0.354824 

0.949248 

0.883022 

0.269637 

0.712632 

0.93301? 

0.176665 

0.450324 

0.969846 

0.091629 

0.216717 

0.992404 

0.031655 

0.056755  Number  of 

Blades  5 

1.000000 

1.010000 

0.000000  Design  Jv 

-  0.889 

Unpublished  experimental  results  at  DTNSRDC  on  a  similar  geometry  indicate  that  the  flow 


over  the  blades  will  be  fully  turbulent  at  a  0.7-radius  Reynolds  number  of  4  \  IOC  A  0.7-radiu 
Reynolds  number  is  defined  as: 

•oo.t  v\/ 1  *  Pr1)2 


R„  = 


(20' 
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FRACTION  OF  RADIUS  XD  FRACTION  OF  RADIUS  X 


V 


where:  (c)q  7  =  blade  chord  at  0.7  radius 

V  =  onset  speed 

Jv  =  advance  coefficient,  Jv  =  V/(nD) 

n  =  constant  rotational  speed,  revolutions 

per  unit  time 
I)  -  rotor  diameter 

y  kinematic  viscosity 

The  assumed  flow  conditions  for  turbulent  flow  are  v  -  1.191  \  l()'6m-/s  (1.2X2  \  10-s  ft-  sec). 
V  =  16.9  m/s  (55.3  ft/sec),  and  a)  =  391.5  rad/s. 

No  previous  analytical  or  experimental  boundary-layer  results  exist  for  these  specific  blades. 
Initial  conditions  are  assumed  to  be  determined  by  flat -plate  values  for  turbulent  flowl-N,29.30it 
with-  0,,  -  1/  U  0.0142  (Rx)6/7,  Cf  =  0.310/ (l»2(2R0)  +  2ln(2Re)],  and  H  - 
I  1 1  -4.67(C'(') 1  -],  and  by  the  rotating  disk  limiting  streamline  angle  ft  value*4*,  ft  -  *  1 1 

degrees  on  the  blade  suction  surface  and  fi  -  -  1 1  degrees  on  the  blade  pressure  surface. 
Because  of  the  uncertainty  in  these  starting  values,  the  computed  results  are  of  qualitative,  but 
not  necessarily  quantitative,  validity. 

The  computed  potential-flow  streamlines  for  the  blades  are  shown  in  figure  S.  The 

streamlines  for  the  warped  Blade  2B  and  the  skewed  Blade  2C  coincide  and  both  are  shown  by 

the  dashed  line  in  Figure  8.  Streamlines  at  increased  radial  values  lie  noticeably  inboard  of  the 

constant  radius  lines  on  the  suction  surface.  F.xcept  near  the  tip,  the  same  trend  is  evident  on  the 

pressure  surface.  Streamlines  13  and  14  coincide  on  the  pressure  surface  of  the  three  blades. 

I  lie  computed  boundary -layer  characteristics  on  the  three  blades  show  similar  trends  along 

each  streamline.  Therefore,  the  results  will  be  presented  only  along  one  streamline,  streamline 
14.  Figures  9  through  II  show  the  computed  distributions  of  8 jj,  tan/T  and  11,  respectively.  In¬ 
creased  values  of  tlie  streatmvi.se  momentum  thickness  Ajj  and  the  shape  parameter  H  are 
predicted  for  the  warped  and  skewed  blades  relative  to  the  unskewed  blade. 

Figure  12  shows  the  predicted  distribution  of  the  local  skin  friction  coefficient  (  t,  as  com¬ 
puted  by  empirical  relationships  given  in  Fquations  (8),  for  Blades  2A,  2B,  and  2(  .  I  he  Mgmti 
cant  result  shown  in  this  figure  is  that  the  local  skin  friction  coefficient  is  predicted  to  be  tcJuc 
ed  on  both  the  suction  and  pressure  surfaces  when  the  blade  is  highly  warped  or  skewed.  Rednc 


-  Rx  (chord)  (fraction  of  chord)  U-v 
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lions  of  up  to  30- percent  are  predicted  near  the  trailing  edge.  Similar  reductions  are  predicted 
for  each  streamline.  Small  variations  in  the  predicted  values  of  the  potential  How  velocity  for  the 
various  blade  designs  do  not  significantly  effect  these  results.  Since  insufficient  data  are  com¬ 
puted  to  permit  a  calculation  of  total  blade  drag,  an  average  reduction  of  1 5-percent  was 
employed  for  the  average  ehordwise  blade-section  drag  relative  to  a  constant  drag  coefficient  of 
0.0085.  Thrust  loading  coefficient  and  power  coefficient  increased  by  one-half-percent  and 
reduced  by  one-percent,  respectively,  hence  an  increase  of  efficiency  of  about  two-percent  for  a 
warped  or  skewed  blade.  The  actual  change  in  efficiency  was  1  percentage  point. 

Several  calculations  were  made  to  identify  the  factors  contributing  to  the  predicted  reduc¬ 
tion  in  skin  friction  with  increased  blade  warp  or  skew.  C  alculations  to  date  seem  to  indicate 
that  these  differences  arise  from  geometrical  effects.  Relative  to  the  unskewed  blade,  the  poten¬ 
tial  flow  velocity  decreases  by  less  than  2  percent  and  the  metric  coefficient  hj  decreases  by  10  to 
20  percent  for  both  the  highly-warned  and  the  highly-skewed  blades.  The  metric  coefficient  h-» 
inc; eases  bv  10  to  15  percent  for  the  warped  blade  and  by  2  to  3  percent  for  the  skewed  blade. 

A  reduction  is  predicted  in  the  angular  velocity  component  of  S)  to  14  percent  for  the  warped 
blade  and  of  only  l  to  3  percent  for  the  skewed  blade.  Unis,  ii  does  not  appear  as  though  ativ 
single  factor  is  responsible  f*»,  the  predicted  skin  friction  reduction  with  increased  blade  warp  or 
'kew .  Instead,  the  interaction  of  various  geometrical  factors  seems  to  contribute  to  this 
i  eduction . 

<  ( >\C  1  ID! NC i  R1  MARKS 

1  his  paper  presents  analysis  and  results  for  computing  three-dimensional  turbulent  houn 
darv  layers  on  the  surface  of  a  rotating  blade  using  a  momentum-integral  method.  The  analytical 
and  numerical  procedures  for  solving  the  momentum-integral  boundary-layer  equations  have 
been  presented  and  implemented  in  an  efficient  computer  code.  Typical  boundary-layer  com 
putaiions  require  less  than  It)  seconds  C  PU  on  the  high-speed  Burroughs  7700  computer  ai  a 
cost  less  than  one  dollar.  The  predicted  boundary-layer  parameters  are  shown  to  give  reasonable 
agreement  with  experimental  data  for  a  simple  body.  Assuming  initial  boundary-layer 
characteristics  determined  by  flat-plate  values  for  turbulent  flow,  increasing  blade  uaip  is 
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predicted  to  decrease  the  local  skin  friction  coefficient.  I  he  reduction  in  skin  friction  increases 
the  thrust  coefficient  and  reduces  the  torque  coefficient  yielding  an  increased  elfieiencv  ol  about 
two  percent  (one  percentage  point).  However,  the  accuracy  of  the  initial  \ allies  and  the  effect  ol 
their  variation  on  the  predicted  results  is  not  known. 

Future  work  in  this  area  should  include  an  improved  crossflow  veloeitv  profile,  belter 
estimation  of  the  initial  conditions  for  the  turbulent  boundarv -layer  calculation,  and  the  caiculu 
tion  of  thrust  and  torque  to  predict  drag.  Fxperimental  boundarv -layet  data  are  needed  to  inllv 
assess  and  improve  the  present  analytical  work.  Future  work  should  also  include  the  prediction 
of  laminar-flow  characteristics  at  the  leading  edge,  laminar  to-turbulence  ttansistion  piedictioiis. 
and  eventually  the  prediction  of  separation  phenomena  both  on  the  blade  uilaee  and  at  the  re 
where  a  lip-vortex  forms.  Pioneering  work  on  the  tip  tcgion  lias  been  published  hv  slamuoih.  e: 
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The  turbulent  boundary  layer  results  for  Blades  2  (unskewed  Blade  2A,  72  degree 
warped  Blade  2B,  and  72  degree  skewed  Blade  2C)  presented  in  the  main  text  are  com¬ 
puted  for  a  model-scale  0.7  radius  Reynolds  number  of  4  *  10^.  In  this  appendix, 

similar  results  are  presented  for  these  blades  for  a  full-scale  0,7  radius  Reynolds 

7  -6  2 

number  of  6  *  10  .  The  assumed  flow  conditions  are  V  =  1.191  *  10  m V s  (1.282  x 

-5  2 

10  ft^/sec),  V  =  16.9  m/s  (55.3  ft/sec),  and  u)  =  26.1  rad/sec.  The  propeller 
diameter  is  4.572  m  (15  ft). 

A  finer  computation  grid  was  selected  for  the  full-scale  computations.  More 
grid  nodes  were  added  in  the  radial  direction  in  order  to  compute  the  boundary  layer 
characteristics  closer  to  the  hub  at  the  trailing  edge  of  the  blade.  The  boundary 
layer  results  on  the  model-scale  blades  were  also  recomputed  using  the  finer  grid. 
Both  full-scale  and  model-scale  results  are  given  in  this  appendix. 

Again,  initial  conditions  for  the  streamwise  momentum  thickness  u  and  the 

sliape  factor  H  are  assumed  to  be  determined  by  flat-plate  values  for  turbulent 
28,29  30 

flow  as  described  in  the  main  text.  The  value  of  the  limiting  streamline 

angle  £  is  also  given  by  the  rotating  disk  value  of  ft  =  +11  degrees  on  the  blade 
suction  surface  and  £  =  -11  degrees  on  the  blade  pressure  surface.  Since  no  previ¬ 
ous  analytical  or  experimental  data  exist  for  these  blades  and  the  starting  condi¬ 
tions  are  at  best  a  guess,  the  computed  results  can  be  considered  only  qualitatively 
val id . 

Since  the  local  skin  friction  coefficient  C ^  contributes  to  the  overall  drag, 
it  is,  perhaps,  the  boundary  layer  result  of  most  interest.  Therefore,  the  only 
comparison  presented  is  that  of  the  skin  friction  coefficient  at  the  0.73-radius  for 
the  three  blades  at  both  model-scale  and  full-scale.  Figure  13  shows  the  computed 
percentage  change  in  C^.  for  Blades  2B  and  2C  from  the  unskewed  Blade  2A.  The  open 
symbols  represent  the  full-scale  results  and  the  solid  symbols  represent  the  model- 
scale  results.  As  shown  in  Figure  13a,  a  decrease  in  the  skin  friction  coefficient 
is  predicted  on  the  suction  surface  at  both  scales  for  the  warped  and  skewed  blades. 
The  decrease  in  continues  to  the  trailing  edge  of  the  blades.  There  is  also 
1  Lttle  variation  in  the  results  for  warp  and  skew  at  each  scale.  Slightly  larger 
changes  in  C^.  are  predicted  for  the  model -scale  blades  than  for  the  full-scale  blade 
Figure  13b  shows  the  changes  in  the  local  skin  friction  coefficient  on  the  pressure* 
surface  of  the  blades.  Again,  a  decrease  in  C^.  Is  predicted  at  both  full-scale  and 
model-scale  for  blade  warp  and  blade  skew.  However,  at  the  trailing  edge*  of  the 


f 


blade,  varies  little  from  the  value  predicted  for  the  unskewed  Blade  2k.  Blade 
skew  is  shown  in  Figure  13b  to  have  a  smaller  effect  on  C^.  than  blade  warp.  Overall, 
Figure  13  shows  that  blade  warp  and  blade  skew  are  predicted  to  decrease  the  local 
skin  friction  coefficient  C^ . 
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